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ABSTRACT 



Context. Nonthermal radio emission in massive stars is expected to arise in wind-wind collisions occurring inside a binary system. 
One such case, the O-type star CygOB2#9, was proven to be a binary only four years ago, but the orbital parameters remained 
uncertain. The periastron passage of 201 1 was the first one to be observable under good conditions since the discovery of binarity. 
Aims. In this context, we have organized a large monitoring campaign to refine the orbital solution and to study the wind-wind 
collision. 

Methods. This paper presents the analysis of optical spectroscopic data, as well as of a dedicated X-ray monitoring performed with 
Swift and XMM-Newton. 

Results. In light of our refined orbital solution, Cyg OB2 #9 appears as a massive O+O binary with a long period and high eccentricity; 
its components (05-5.51 for the primary and 03-4111 for the secondary) have similar masses and similar luminosities. The new data 
also provide the first evidence that a wind-wind collision is present in the system. In the optical domain, the broad Ha line varies, 
displaying enhanced absorption and emission components at periastron. X-ray observations yield the unambiguous signature of an 
adiabatic collision because, as the stars approach periastron, the X-ray luminosity closely follows the 1/D variation expected in that 
case. The X-ray spectrum appears, however, slightly softer at periastron, which is probably related to winds colliding at slightly lower 
speeds at that time. 

Conclusions. It is the first time that such a variation has been detected in O+O systems, and the first case where the wind-wind 
collision is found to remain adiabatic even at periastron passage. 

Key words, binaries: spectroscopic - stars: early-type - stars: emission-line - stars: individual (Cyg OB2#9) - X-rays: stars 



1. Introduction 

Nonthermal radio emission, bright and variable hard X-ray emis- 
sion, or periodic changes in the profiles of optical emission lines 
are the typical signatures of colliding winds in massive binary 
systems. Indeed, only massive stars have strong enough stellar 
winds for hosting these kinds of energetic phenomena. However, 
the reverse is not true; i.e., (many) massive binaries exist that 
show no signature of colliding winds. In fact, nonthermal radio 
emission is a rare feature (fewer than 40 cases in our Galaxy, see 
iDe Becke3l2007l; l Be naglia 2010), and recent X-ray surveys have 
shown that most bi naries are not significantly overluminous at 
high energies (Naze 2009). Understanding why it is so will im- 
prove our knowledge of both collision mechanisms and stellar 
wind properties since such emissions heavily depend on them. 
To this aim, it is important to collect as much information as 
possible, notably on the rarest cases showing conspicuous colli- 
sion signatures like Cyg OB2 #9. 

As suggested by its name, this star belongs to the Cyg OB2 
association, one of the richest OB associati ons of our Galaxy 
since it contains about one hundred O-stars (Knodlseder 2000; 
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Negueruela et al. 2008; Wright et"aTll2010t) . As for other mem- 
bers of this association, CygOB2#9 is strongly extinguished, 
which renders its study q uite difficult (E[B - V] = 2.11, 
Massey & Thompson 1991). It was nevertheless soon identified 
as an early-type star (Q5If + . lMorgan et alJIl 955b IWalbornll 19731) 
and as a nonthermal radio emitter dAbbott et al.l I1984TT After 
these first studies, the spectrum of Cyg OB2 #9 was not studied 
in more detail, so that it was only proven to be a binary in re- 
cent years. A hint towards a multiple nature came from the conic 
shape of the nonthermal emissio n, a configuration reveale d by 
the Very Large Baseline Array (Dougherty & Pittard 2006) and 
only possible in wind-wind collisions. Two years later, two im- 
portant findings were reported. A first, indirect evidence of bina- 
rity came from the periodic modulations of the radio emission 
(Ivan Loo et al.l 120081) . occurring with a timescale of 2.355 yr. 
Soon after that, direct evidence of binarity, i.e., clear doubling 
of stel lar lines, was unv eiled by a dedicated spectroscopic cam- 
paign dNaze et al.l l2008). The continued monitoring of the sys- 
tem then led to c omputation of a preliminary orbital solution 
dNaze et al.ll2.010h : CygOB2#9 appeared as a long-period, very 
eccentric (e = 0.7) binary composed of two nearly equal-mass 
objects. In parallel, the analysis of several XMM-Newton ob- 
servations of this object revealed no strong overluminosity or 
peculiar hardness in the X-ray spectrum, but did show some un- 
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usual long-term variability (by about 10%, iNaze et aT] 12010). 
Clearly, additional data are needed to pinpoint the properties of 
the system: improving the orbital solution, enlarging the orbital 
coverage of the X-ray data (the available ones being all taken 
far from periastron), and extending the knowledge of the radio 
emission. Only then can a full modeling be undertaken, to better 
understand massive stars and their winds. 

Since the discovery of the binarity of Cyg OB2 #9, two pe- 
riastron passages have occurred. The first one took place in the 
beginning of 2009, when the system was in conjunction with 
the Sun and therefore unobservable. The second one was ex- 
pected to occur in June-July 2011, i.e., at the best time for ob- 
serving the Cygnus constellation. We thus organized a large ob- 
serving campaign of Cyg OB2 #9 to observe all aspects of this 
event: radio monitoring of the nonthermal emission using the 
Expanded Very Large Array, interferometric monitoring to dis- 
entangle the astrometric orbit of the two stars using CHARA, 
spectroscopic monitoring in the optical to improve the orbital so- 
lution, and X-ray monitoring to detect and understand the high- 
energy emission from the wind-wind collision. The results of 
this multiwavelength campaign will be published as a series of 
papers, which will be concluded by a modeling of the system. 
With all these observations in hand, CygOB2#9 will become 
one of the best known nonthermal radio emitters, with detailed 
constraints on the wind characteristics and on the collision prop- 
erties. CygOB2#9 will thus represent a true stepping stone in 
our understanding of wind-wind collisions and their associated 
emissions. 

This first paper provides the results obtained through optical 
spectroscopy and X-ray monitoring. Section 2 presents the ob- 
servations, Sect. 3 analyzes the spectroscopic data with the aim 
to improve the orbital solution, Sect. 4 discusses the variations 
observed in X-rays, and Sect. 5 summarizes our findings and 
concludes this paper. 

2. Observations 

2. 1 . The optical domain 
2.1.1. Sophie 

Our monitoring of CygOB2#9 at the Haute-Provence 
Observatory (France) began in 2007. It mainly used the 
high-resolution echelle spectrograph Sophie installed on the 
1.93m telescope. It yields an homogeneous dataset covering a 
broad wavelength domain with the highest quality available. 

The Sophie spectrograph provides coverage of the 3900- 
6900 A domain with a resolving power R = A/AA of 35 000. 
In total, 75 echelle spectra were obtained with this instrument. 
The exposure time was typically 40 min, and the derived signal- 
to-noise ratios vary from 30 to 80. Whenever the exposures were 
close (i.e., they were taken within a month and no large radial- 
velocity [RV] variations were detected), these data were com- 
bined, providing in the end 31 independent spectra spanning a 
period of four years. As for all Sophie spectra, these data were 
reduced by the automatic Sophie pipeline Q it locates the or- 
ders, performs an optimal extraction with cosmic -ray rejection, 
calibrates the spectra, and finally reconnects the orders after 
barycentric correction and correction for the blaze function (de- 
rived from a Tungsten lamp exposure - for massive stars, it pro- 
vides only a rough, first-order normalization). These data were 
smoothed further using a sliding window of 15 px size (with 

1 http://www.obs-hp.fr/guide/sophie/sophie-eng.shtml 



1 px=0.01 A) , and corrected by hand for any remaining cosmic 
ray hits. The journal of observations is provided in Table Q] 

Owing to the low signal-to-noise and the high extinction of 
the star, it is difficult to perform a clean and coherent normal- 
ization. To ensure homogeneity, we first performed the best nor- 
malization possible, using low-order polynomials, on the data 
from June 2010: these data are not only among the best qual- 
ity data, but they were also taken at a time when lines are most 
blended, making the definition of continuum windows easier. We 
then used this normalized set as a reference for the other data. It 
enables us to remove stellar lines (through a division of the spec- 
tra by the reference), allowing better modeling with low-order 
polynomials of the continuum oscillations due to the poor blaze 
correction and poor order merging. 

In a final step, we tried to correct for some annoying non- 
stellar features. To this aim, we have fitted with Gaussians the 
diffuse interstellar bands (DIBs) in the 5412 A and 5800 A re- 
gions (see Fig. UJ on the June 2010 data (i.e., when lines are most 
blended). In addition, a DIB near 58 12 A was unveiled when the 
stellar lines were deblended - this particular DIB apparently has 
a similar rest wavelength to the C iv line. It was thus fitted on 
the deblended spectra and added to the set of DIBs already de- 
fined. No good fitting could be found for the DIB at 5797 A, so 
that we kept it in the spectra. The fitted interstellar features were 
then removed from the data, l eading to cleaner st ellar lines. We 
have also used the table from lHinkle et alj d2000h to correct the 
Ha, Heil5876A, and Cml5696A regions for pollution by the 
telluric lines. 

2.1.2. WIRO 

Our team also monitored Cyg OB2 #9 at the Wyoming Infrared 
Observatory (WIRO) 2.3m telescope many times between 201 1 
May 31 and 2011 July 25 (see Table Q]), which provides a 
more intense coverage of the periastron event, as well as an 
independent check of the stellar behavior at that time. All 
new datasets were reduced using standar d IRAF reduction rou- 
tines as outlined in iKiminki et al.l d2007l) . Sixteen observations 
were obtained using the WIRO-longslit spectrograph with the 
1800 line mirT 1 grating in first order over the spectral range 
5250-6750 A. Exposure times varied from 720 s to 5400 s de- 
pending on weather conditions; the resulting signal-to-noise was 
about 200 near 5876 A. Copper- Argon lamp exposures were 
taken after each observation to wavelength-calibrate the spectra 
to an RMS of 0.03 A. The typical resolution was 1 .5 A FWHM 
across the chip. Spectra were Doppler corrected to the heliocen- 
tric frame and then shifted further by small amounts (typically 
< 6kms _1 ) so that Nai interstellar absorption lines in each ex- 
posure fell at the mean Na i velocity averaged over the ensem- 
ble of 16 measurements. The shift also ensures that interstel- 
lar features, e.g. the Nail 5889 line, display similar velocities 
in WIRO and Sophie spectra, enabling us to directly compare 
them. Finally, continuum normalization was performed using 
low-degree polynomials. 

2.2. X-rays 
2.2.1. Swift 

The Swift facility observed Cyg OB2 #9 during 5 ks at five care- 
fully selected dates in 2011 (PI Naze, see Table [TJ, to follow 
the star during ten months around the date of periastron. The 
data were then reduced following the recommendations of the 
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UK Swift centefl First, the pipeline (FTOOLS task xrtpipeline) 
was applied. Events were then extracted in a circular region of 
47" radius for the source and in a nearby circle of 119" radius 
for the background. T he chosen re gion is as devoid of X-ray 
sources as possible (see Rauw j20i"ll) . The adequate response ma- 
trix file (RMF, swxpc0tol2s6^0010101v013.rmf) is provided 
by the Swift team, but the ancillary response file (ARF) was 
specifically calculated for Cyg OB2 #9 by the task xrtmkarf with 
the inclusion of an exposure map so that bad columns can be 
taken into account. Count rates measured in these observations 
vary from 0.04 to 0.1 cts s , so there is no pile-up. 

2.2.2. XMM-Newton 

XMM-Newton had observed the Cygnus OB2 region six times 
in previous years, mostly for studying the X-ray emission of 
Cyg OB2 #8A (see TableQ]). Since Cyg OB2 #9 was in the field- 
of-view, we analyzed its X-ray emissi on and found some un- 
explained variations dNaze et al.l 1201 Oh . A full monitoring of 
CygOB2#9 with XMM-Newton in 2011 was not possible, but 
an observation was scheduled near the predicted time of the pe- 
riastron passage (25 ks, ObsID=0677980601, PI Naze, see also 
TableQ]). 

This new dataset was reduced in the standard way using SAS 
v. 10.0.0. After the pipeline processing (SAS tasks emproc, ep- 
proc), we filtered the data to only keep the best events, as rec- 
ommended by the SAS teanfl A flare in the X-ray background 
affects the middle of the observation. It was eliminated, so that 
the effective exposure times are reduced to 24 ks for MOS and 
15 ks for pn. 

Source events were extracted in a circular region of 35" ra- 
dius around the target, and the background region was defined 
as a nearby circle of the same size. An annulus could not be 
used, as num erous X-ray sources can be found in the region (see 
lRauwll20Tlh . The count rate measured for this new observation 
is 0.6 cts s for MOS and 1.8 cts s _1 for pn, which is five times 
higher than in previous observations. Despite this brightening, 
the source is not affected by pile-up, as checked by a run of the 
task epatplot. 

3. Cyg OB2#9, as seen in the optical 

3. 1 . The spectra 

The high extinction of Cyg OB2 #9 has two direct consequences 
on its optical spectrum. First, the blue part of the spectrum 
(< 4400 A) is very noisy, hence not really usable. Second, most 
lines recorded in the spectra actually are interstellar features! 
Only a few stellar lines can be distinguished and used for further 
study (see Fig.©: Balmer lines (Ha, HyS), Heil4471, 5876A, 
Hen 14542, 4686, 5412, 6683 A, N ml 4634, 4641 A (in emis- 
sion),^ ra 1 4650, 5696, 6723, 6729A (in emission), Civ 15801, 
5812A. Some of these stellar lines are, however, contaminated 
by interstellar features, such as Civl5812A. 

When one compares the spectra taken at different phases (see 
Fig. [3]), it is obvious that the Cral5696A emission line follows 
the motion of one component of the system. The He n 1 4686A 
emission line instead displays the opposite motion, that of the 

2 http://w w w. swift, ac.uk/analysis/xrt/ 

3 See threads available on 
http://xmm.esac.esa.int/sas/current/documentation/threads We kept 
single and double events for MOS (PATTERN=0-12 <m&XMMEA_EM 
filter) and single events for pn (null PATTERN and FLAG=0-4). 
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Fig. 3. Comparison of the spectral appearance of CygOB2#9 
in June 2010 (0=0.571, solid line) and in the beginning of July 
2011 (0=0.005, thick dotted line). The spectra have been cor- 
rected for some DIBs (see text) and the telluric lines (the latter 
correction being imperfect around 5696 A). 

other component, though with a shift of about 90kms~'. This 
shift implies that the line most probably does not arise from the 
photosphere, but somewhere close to it, in the wind. The exact 
position may be influenced by the presence of a wind-wind col- 
lision in the system (see below). These two lines are the sole 
ones "associated" to a single component of the system. Other 
lines show two components in the 0.9-1.1 phase interval (Fig. 
[3}. Another problem is obvious when one looks at the spectra 
where both components are clearly separated: the lines of both 
objects have similar strengths in He n, N m, and C iv, which does 
not help the derivation of the RVs and their attribution to one of 
the components. 

3.2. RV measurements and orbital solution 

As already explained in Naze et al. (2010), measuring the RVs 
in CygOB2#9 is difficult because the broad stellar lines are 
faint, sometimes even contaminated by interstellar lines, and 
also blended for most of the (highly eccentric) orbit. In a first 
step, we measured the RVs by fitting one (or two, if the line dis- 
played two components or appeared asymmetrical, which hap- 
pens within A<p = 0.1 of periastron) Gaussian(s) to the stellar 
lines. Results of these simple fits are shown at the top of Fig. [4] It 
is obvious from that figure that deblending the two components 
is not easy, since they can only be measured at maximum sepa- 
ration. The shift of He n 1 4686A with respect to the true stellar 
motion is also obvious, as is the attribution to the primary star of 
15696AQ and the main component of He i lines. WIRO data 

4 The WIRO measures of this line around 0=0. 9 are very uncertain 
due to the unexplained presence of a faint, broad second component 
(normalization problem?). 
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Fig. 2. Spectrum observed in June 2010 (after DIB and telluric correction). The labels indicate the stellar lines (other features are 
interstellar). The He n 1 6560A line next to Ha is strongly blended with a DIB. 



agree well with the Sophie observations, when taking their larger 
uncertainty (about 20-30 km s _1 ) due to a lower resolution into 
account. 

To improve over these simple fits, we fixed the component's 
shapes to the Gaussians fitted in the most separated cases, and 
we cross-correlated the m with the b l ended lines until the best^- 2 
was found, as we did in Naz e et al. (2010). The new RV curves 
are better, as can be seen in Fig. [4] However, a trend must be 
noted: the x 1 mapping has a tendency to separate the compo- 
nents more at phases of maximum blending. This may be due to 
slight changes in line shapes (due to noise and imperfect normal- 
ization), but allowing the line shapes to vary slightly produces 
erratic results of lower quality. 

Finally, we used these preliminary RV s as input for a dis- 
entang ling program based on the method of Gonzalez & Levato 
(2006). In this program, the spectra of each component are it- 
eratively determined, first on the basis of the approximate input 
RVs, then by using improved RVs found from cross-correlating 



the components' spec tra with the observed spectra (for details 
see iMahv et all 1201 Oh . In our case, the disentangling tries to 
fit the spectra of three components: the fixed interstellar fea- 
tures, and the primary and secondary stars (whose RVs vary 
from one observation to the next). The continuum was arbi- 
trarily fixed at 0, 0.5, and 0.5 for these three components (i.e., 
the true light ratios are not taken into account in this method). 
The cross-correlation was performed using masks, containing 
the Heil4471A, Henl4542A, Henl5412A, and Heil5876A 
lines for the primary and the Hen 14542 A, Hen 154 12 A, and 
Heil5876A lines for the secondary. The spectra where both 
components are clearly visible (i.e., maximum separation) re- 
ceived a weight 50 times greater than in the blended cases. This 
process yields the RV curves seen in the last panel of Fig.|4](for 
values, see Table 2), and the disentangled spectra are shown in 
Fig.0 

The RVs derived from the disentangling can be used to find 
the orbital solution. To this aim, we used the Liege Orbital 
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Table 3. Orbital solution for CygOB2#9. T (hence , 
responds to periastron passage, in HJD - 2 450 000. 



0) cor- 



Parameter 


Value 


P(d) - from radio 


858.4 ±1.5 


To 


4023.821 ±2.51 


e 


0.713 ±0.016 




192.1 ±2.9 


Mi /M 2 


1.13 ± 0.08 


/l , — 1 \ 

yi(kms ) 


1 1 O i 1 H 

-33.S ± 2.7 


7 2 (kms-') 


0.9 ±2.8 




60.3 ± 3.0 




68.0 ± 3.4 


a i sirwXRo) 


716.9 ±39.5 


a 2 sin /(Rq) 


808.4 ± 44.5 


Mi sin 3 ('(Mo) 


34.3 ± 4.6 


M 2 sin 3 ((Ms) 


30.4 ± 4.0 


RMS(kms-') 


9.6 



Solution Package ( LOSPfl which is based on the SB1 method of 
I Wolfe etail (|1967|) g eneralized to SB2 cases as described, e.g., 
in lRauw et al.l(l2000al) . The period was fixed to 858. 4d, as found 
from an extensive radio dataset covering three decades (Blomme 
et al., in preparation^ and the derived parameters are shown in 
Tableland Fig(6]. It agrees well with the preliminary solutions 
from lNaze et al.fl2010l) . 

3.3. Spectral features 

While photospheric Hen absorption lines display similar 
strength in both components (Fig. |3J, the Hei lines appear 
much stronger in the primary component (i.e., the one that 
appears blueshifted at periastron). In fact, Heil4471A is so 
weak in the other component that it is difficult to distin- 
guish it. Comparing the relative strength of Hen 14542 A and 
He 1 14471 A dContill 19731) . we found that they correspond to an 



05-5.5 typ e and Q3-4 type for the primary and secondary, re- 
spectively dWalborn et alj|2010l reported the blended spectrum 
as 04.5Ifc). That one component is hotter than the other one 
may explain why C m 1 5696 is best seen in the cooler compo- 
nent while C iv lines have similar strengths in both stars. Since 
the He n lines display similar intensities, it also suggests that the 
two stars should have different luminosity classes, the coldest 
object being also bigger. This conclusion is supported by the 
fact that identical luminosity classes would also be incompati- 
ble with the coldest star being the most massive component (see 
previous section). 

To investigate this issue further, we turned to other spec- 
tral features. Both components display Nm and the neighbor- 
ing C m lines in emission. The latter lines are not variable and 
are much stronger in the primary, which agrees well with the 
finding th at these features r each their maximum strength in 
05 stars dWalborn et al.ll2010h . In our disentangled spectra, the 
He ii 14686 A line appears in emission in both components, but 
the width of the feature in the primary spectrum suggests that 
this line may be contaminated. Indeed, the disentangling pro- 
gram is not fit to study such features since it considers the line 
shapes to be constant throughout the orbit and to arise solely in 
the stellar atmospheres. The shift observed for the Henl4686A 



5 The LOSP package is maintained by H. Sana. It is available, with a 
preprint describing it, from http://staff.science.uva.nl/~hsana/losp.html 

6 The time interval covered by the spectroscopic data, about 1500 d, 
is too short for a precise period determination. 



line (see above) clearly indicates that it is not the case: even if the 
stars might truly display some He n 1 4686A in emission, most of 
it should be formed in the wind, close to the secondary star and 
with a potential influence of the wind-wind collision. With N m 
+ Cm lines (and possibly Hen 14686 A) in strong emission, the 
primary spectrum thus appears quite compatible with a super- 
giant class of type fc. 

Determining the luminosity class of the secondary star is 
more difficult. We ch ecked the spectra for signs of N iv 14058 A 
dWalborn et al.l l2002h : this line appears as a very weak emis- 
sion, much weaker than N ml 4634,4641 A, but the spectra are 
so noisy in the blue range that it is difficult to perform disentan- 
gling, hence to say whether the line belongs to one component 
or to both. In either case, this suggests that the hotter compo- 
nent (the secondary) is 03.5V or 04III-V and not a supergiant, 
whereas the presence in its spectrum of Nm and Hen 14686 A 
lines in emission suggests the contrary! 

We used the typical s tellar parameters of giant 03.5 and 05.5 
stars dMartins et al.ll2005l temperature and gravity are similar 
for giants and supergiants in this reference) as input for the atmo- 
sphere modeling code CMFGEN. This enabled us to determine 
the equivalent widths of the He 1 1 447 1 A and He n 1 4542, 54 1 2 A 
lines on the synthetic spectra, which we compared to the values 
found by fitting Gaussians to the observed spectra (at maximum 
velocity separation) or to the disentangled spectra. This led to 
nearly equal luminosity ratios, while the orbital solution further 
suggests equal masses. Therefore, our favored spectral types are 
05-5.51 for the primary and 03-4III for the secondary. 

Our monitoring of Cyg OB2 #9 allows us to search for spec- 
tral signatures of a wind-wind collision. To this aim, one gen- 
erally turns to Balmer and Henl4686A emissions. As already 
mentioned, the Henl4686A line has a peculiar behavior: its 
shape is relatively constant, but its RVs are shifted relative to 
the stars (Fig. 3]). The recorded velocities instead suggest the 
line arises close to the secondary star (but not close to its pho- 
tosphere, because of the velocity shift), and a contamination by 
the collision zone is possible but difficult to assess. 

In contrast, the broad and intense Ha emission line displays 
no obvious shift with orbital phase but its profile clearly varies 
though the changes are not dramatic (Fig. [7]). We first describe 
these variations. In this context, one has to keep in mind that 
only a single periastron passage could be observed in details, 
that of 2011, so that caution must be applied when generaliz- 
ing the observations with </> ~ to the usual behavior of the 
star. At phases <p > 0.1, the line displays an overall triangu- 
lar shape with a basis covering about 2000 km s _1 and a peak at 
~200kms~'. Superimposed on this emission are two relatively 
narrow absorptions with velocities of about -175 and 45 km s _1 
(i.e., not stellar velocities). It should be noted that the RVs of 
these two absorptions do not seem to vary much, but neither do 
the stellar RVs at these phases. As time goes by, the emission 
peaks decrease slightly in strength, though, as do the absorptions 
(the blueshifted absorption disappearing only after the redshifted 
one). Before periastron ((/> = 0.85 - 0.98), the line thus appears 
as a nearly perfect triangle. Suddenly, shortly before periastron 
(at about </> = 0.99 to the best of our knowledge - our data do 
not enable us to further pinpoint the exact time of the event), the 
emission peak becomes more intense and is shifted to the red. 
It now appears at about 300kms _1 . In parallel, a strong absorp- 
tion appears at ~ -200 km s -1 , and is also much broader than 
before, since it is approximately 60% wider than the previous 
blueshifted absorptions. A substructure of this absorption may 
be related to the photospheric absorption by the primary star, 
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though this component is not obviously detected at other phases. 
In the WIRO data, the emission appears strongest at (f> — 0.999, 
whereas the absorption is most intense at <p — 0.998, but this can- 
not be confirmed by Sophie data owing to the lack of coverage 
of these phases. After periastron, the emission peak decreases in 
strength and moves to less blueshifted velocities, while the ab- 
sorption becomes less broad, and a second absorption gradually 
appears, coming back to the start of the cycle. 

A very broad Ha emission is quite compatible with a line 
arising not at the photosphere, but over a wide area - maybe 
the wind-wind collision in our case. If the wind-wind collision 
zone is dense enough at periastron, it could absorb some light 
from the stellar photospheres, in a kind of atmospheric eclipse. 
In such a case, the broad absorption appearing at periastron may 
be related to absorption of the primary's light by the colliding 
wind zone (Fig.[6]l. Indeed, the collision zone ma y well be h ighly 
distorted at periastron by Coriolis effects (e.g. iParkin & Pittardl 
2008), and one of its "wings" could then obscure our view of the 
primary, as seen from Earth. In such a configuration, the emis- 
sion would come from the other zones of the collision, which 
are mostly redshifted and do not appear in front of one star, as 
seen from the observer's point of view. The pair of narrower ab- 
sorptions seen at other phases is more difficult to explain, since 
the density of the collision zone decreases strongly with increas- 
ing separation, but they could be explained by a similar effect if 
there is enough dense cold gas downstream. Their constant ve- 
locity, though, is not surprising since the stars themselves do not 
show much RV change at these phases. The disappearance of all 
structures at <p — 0.85 - 0.98 is, however, puzzling. It could be 
related to a physical change in the collision zone, or be simply 
due to viewing angle effects. 

4. Cyg OB2#9, as seen in X-rays 

4. 1 . The spectra 

First, we fitted the new XMM-Newton spectra and the old ones 
with optically-thin thermal plasma models of the form wabs * 
phabs * (apec + apec). The fitting was performed within Xspec 
vl2.4.0ad (which uses APEC 1.3.1): the first absorption was 
fixed to the interstellar one, 1.15 x 10 22 cm~ 2 , and the second 
one allows for additional, in-situ absorption by the windfl; two 
thermal emission components were necessary to obtain a good 
fit. Because these XMM-Newton data display the higher signal- 
to-noise ratios, this fitting enables us to see whether temperatures 
and absorption change with phase. The values found for the ad- 
ditional absorption and the cooler temperature in the different 
datasets were found to remain within one sigma of each other, so 
that we decided to fix them. Swift data have fewer counts, hence 
are noisier. A fully free 2T fitting results in large error bars and 
sometimes erratic variations. It was thus not possible to check 
the absorption and cooler temperature constancy on these data. 
Therefore, all Swift spectra were fitted with the cooler temper- 
ature and additional absorption fixed. This seems a reasonable 
assumption since (1) the cooler component has a temperature 
typical of the intrinsic emission of hot, massive stars (see e.g. 
Naze 2009, and references therein), and (2) the additional ab- 
sorption is probably associated to the stellar winds, and their 
observed densities should not change much with phase since 
the two stellar components are similar (unlike what happens in 

7 It may be noted that due to the high interstellar extinction, there 
is no usable data at the lowest energies: using an ionized or a neutral 
absorption for the wind therefore does not change the results, since they 
mostly differ below 1 keV. 



WR+O systems). The contribution from a wind-wind collision 
should occur over a wider temperature range, and thus be the 
sole one responsible for the hotter component, and it should also 
be more variable than intrinsic emissions, especially in highly 
eccentric systems (see below). Table|4]yields the results of these 
fits, and Fig. [8] illustrates the variations of the parameters with 
orbital phase. 

We also tried to fix the hotter temperature further, but the 
X 1 was significantly larger for the XMM-Newton data (Table©. 
However, the flux changes were similar for both XMM-Newton 
and Swift, so that the conclusions reported in the next section are 
not affected by this choice. 

Two datasets were taken close to periastron, 5.5 days apart, 
one by XMM-Newton and one by Swift. They show a differ- 
ence in flux of about 20%. The cross-calibration between Swift 
and XMM-Newton data is believed to be bett er than 20% (A 
Beard more, private communication, see also iTsuiimot o et al.l 
120111) . To check our results, we looked at the other sources in 
the Swift field-of-view. CygOB2#5 and 12 are known to vary, 
but unfortu nately, their variations are not f ully explained and 
predictable dLinder et al.l 12009c lRauwll201ll) . The behavior of 
the colliding-wind binary Cyg OB2 #8A is, h owever, well known 
dDe Becker et al.ll2006t iBlomme et al.ll2010l) , and the Swift data 
are here fully consistent with the XMM-Newton and ASCA data, 
when taking the errors due to a lower number of counts in the 
Swift data into account. The rapid brightening of Cyg OB2 #9 
near periastron, on timescales of days, may thus well be real. 

4.2. Archival data 

The Cyg OB2 association has been observed by all major X-ray 
facilities since its main members were detected in X-rays in 1978 
by Einstein. To see whether our conclusions on Cyg OB2#9 are 
valid, we also considered these old datasets. 

Einstein IPC coun t rates of Cyg OB2 have been reported by 
Wald ron et al.l (I1998L see also the Einstein IPC catalog avail- 
able at HEASARC): they amount to 0.0257 + 0.0018 cts s" 1 
in mid-December 1978 (0 = 0.16), 0.0237 + 0.0042ctss- 1 
on 1979 Nov 22 (</> = 0.56), and 0322 + 0.0014 cts s" 1 on 
1980 Jun 9 {<p = 0.79). However, IChlebowski eTail (fl989l 
HEASARC "ipcostars" catalog) quote an IPC count rate of 
0.0406 + 0.0010 cts s _1 , which is not compatible with the pre- 
vious values though they come from the same data. It is thus 
difficult to know exactly what happened 30 years ago. 

ROSAT observed CygOB2 in April 1991 and 1993 with 
the PSPC (ObsID=200109 and 900314, respectively) and in 
November 1994 and April 1995 with the HRI (ObsID = 201845 
parts 1 and 2, respectively). Spectra of the source (reduced to 
a single bin for HRI) were extracted using the same regions as 
for Swift observations. In these regions, the ROSAT background- 
corrected count rates of Cyg OB2#9 amount to (3.96 + 0.36) x 
10" 2 cts s- 1 in PSPC/1991, to (3.83 + 0.15) x lO^ctss" 1 in 
PSPC/1993, to (2.3 + 0.2) x lO^ctss" 1 in HRI/1994, and to 
(1.47 + 0.06) x lO^ctss 1 in HRI/1995. The P SPC values are 
compatible with the PSPC count rates reported bv lWaldron et al.l 
( 19980 We furthermore fitted the PSPC spectra using the same 
model as above, with temperatures and absorption fixed (the 



8 However, Wald ron et all {T998) report an observed flux in the 0. 1- 
3.5 keV range of ~ 6 x 1CT 13 ergcirT 2 s _1 . We instead found a flux of 
~ 9 x 10~ 13 ergcirT 2 s _I with our best-fit model. The discrepancy may 
come from a typo, as it is well known that ROSAT was not sensitive be- 
yond 2.5 keV: if we take 2.0 or 2.5 keV as the upper limit of the energy 
range, rather than 3.5 keV, Waldron 's value and ours indeed agree. 
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poor signal-to-noise and low sensitivity at high energies did not 
permit us to free more parameters). Fitting results are reported 
in Tableg] 
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Fig. 6. Top: Variation with phase of the RVs (solid line=primary, 
dashed line=secondary), separation, and position angle (defined 
to be zero when the primary is in front and along the line-of- 
sight). Bottom: The orbit of the primary star around the center 
of mass (see Table [3] for physical parameters). The position of 
periastron is shown with a big dot, of phases 0. 1 to 0.9 with open 
squares, of phases 0.05 and 0.95 with crosses, and of phase 0.99 
with a star. The line of nodes corresponds to the x-axis and the 
line-of-sight to the y-axis, with an Earthly observer towards the 
bottom. At periastron, a (arbitrary) sketch of the collision zone is 
shown as a pink dashed line, with the position of the secondary 
star shown as an open blue circle. 
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Fig. 7. Evolution of the Ha line in the Sophie (top) and WIRO 
(bottom) datasets. Phases are quoted near each spectrum and, in 
the top panel, the distance between spectra is proportional to the 
phase difference. Note the two absorptions that disappear before 
periastron, the line profile then being almost triangular, and the 
increased absorption and emission components at periastron. 
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Table 4. Results of the X-ray spectral fitting. 



Date 


HJD 





D/a 




Absorptions and one temperature fixed 






Absorptions 


and temperatures fixed 






-2400000 








normi 


norm 2 


C'obs 
X 


cunabs 
X 


X 2 (dof) 


normi 


norm 2 


pobs 
r X 


rrunabs 


X 2 (dof) 










keV 


f 0- 3 cm" 5 


lO- 3 cm" 5 


lQ- n er£ 


cm 2 s-' 




lO" 3 cm" 5 


lO- 3 cm -s 


10 12 erg cm 2 s 1 




XMM-Newton 




























Oct 04 


53308.580 


0.167 


1.13 


2.57 ± 0.09 


3.56±0.f2 


1.90 + 0.07 


1.75 + 0.04 


5.31+0.11 


Lf6(394) 


2.80 + 0.10 


2.42 ± 0.06 


1.60 + 0.03 


4.85 + 0.10 


1.37 (395) 


Nov 04a 


53318.558 


0.178 


1.17 


2.37 ± 0.09 


3.21+0.12 


1.96 + 0.08 


1.66 + 0.03 


4.96 + 0.10 


1.19 (422) 


2.60 + 0.10 


2.38 ± 0.06 


1.55+0.03 


4.62 ± 0.09 


1.30 (423) 


Nov 04b 


53328.544 


0.190 


1.21 


2.51 + 0.09 


3.80 + 0.11 


1.91+0.07 


1.77 + 0.04 


5.52 + 0.11 


1.28 (417) 


3.13 + 0.10 


2.39 ± 0.06 


1.63+0.03 


5.13 + 0.10 


1.44 (418) 


Nov 04c 


53338.506 


0.202 


1.25 


2.65 ± 0.25 


3.23 ± 0.22 


1.90 + 0.14 


1.73 + 0.06 


5.03+0.18 


1.15 (125) 


2.42 + 0.18 


2.48 + 0.10 


1.57+0.06 


4.54 + 0.16 


1.39 (126) 


Apr 07 


54220.355 


0.229 


1.34 


2.63 ± 0.20 


2.79 + 0.14 


1.70 + 0.10 


1.53 + 0.05 


4.40 + 0.14 


1.08 (214) 


2.17 + 0.12 


2.15 + 0.07 


1.38+0.04 


4.00 + 0.12 


1.29 (215) 


May 07 


54224.170 


0.233 


1.35 


2.52 + 0.16 


2.82 + 0.17 


1.79 + 0.11 


1.55 + 0.05 


4.48 + 0.15 


0.84 (158) 


2.23+0.14 


2.19 + 0.08 


1.41+0.05 


4.10 + 0.14 


0.98 (159) 


Jun 11 


55738.255 


0.997 


0.29 


1.89 + 0.02 


13.8 + 0.22 


8.83+0.13 


6.36 ± 0.03 


21.1 ±0.1 


1.59 (444) 


13.8 + 0.16 


8.83+0.08 


6.36 ± 0.03 


21.1 ±0.1 


1.59 (445) 


Swift 






























Jan 11 


55571.618 


0.803 


1.24 


1.76 + 0.45 


1.51 ± 1.20 


3.04 ± 0.82 


1.62 + 0.21 


4.14 + 0.54 


0.53 (18) 


1.93+0.73 


2.81+0.40 


1.66 + 0.21 


4.32 + 0.56 


0.51 (19) 


Apr 11 


55655.837 


0.901 


0.82 


2.00 ± 0.48 


1.68+1.50 


4.17 + 0.84 


2.41+0.24 


5.47 ± 0.54 


1.45 (25) 


1.47 + 0.85 


4.38 ± 0.50 


2.36 ± 0.23 


5.18 ±0.51 


1.40 (26) 


May 11 


55700.082 


0.953 


0.51 


3.26 ± 0.72 


9.05 ± 1.10 


3.82 ± 0.65 


4.36 ± 0.30 


13.0 + 0.9 


1.06(47) 


6.25+1.10 


5.69 ± 0.58 


3.71+0.25 


11.1 ±0.8 


1.29 (48) 


Jul 11 


55743.839 


0.004 


0.29 


2.15 + 0.38 


10.9+1.7 


6.53 ± 0.98 


5.24 ± 0.30 


16.6 + 0.9 


0.93 (54) 


9.53 ± 1.30 


7.48 ± 0.68 


5.07 ± 0.29 


15.8 ±0.9 


0.93 (55) 


Oct 11 


55842.169 


0.118 


0.91 


6.16 + 3.60 


6.09 ± 1.00 


1.40 + 0.65 


2.50 + 0.38 


7.87 ± 1.20 


0.55 (17) 


4.31+0.99 


2.59 + 0.50 


1.89 + 0.29 


6.44 ± 0.98 


0.78(18^ 


ROSAT 






























Apr 91 


48368.064 


0.41 


1.67 














2.48 ± 1.10 


2.03 ± 1.20 


0.49 ± 0.05 


3.13 ±0.33 


0.65 (lj| 


Apr 93 


49109.312 


0.28 


1.46 














2.93 ± 0.44 


1.30 + 0.44 


0.45 ± 0.02 


3.21 ±0.14 


0.97 (58) 


ASCA 




























e. 


Apr 93 


49106.958 


0.27 


1.45 














1.30 + 0.31 


1.71+0.18 


1.03+0.06 


2.79 ±0.17 


0.96 (19%) 



Notes. Fitted models are of the form wabsxphabsx(apec+apec). The absorptions are fixed to 1.15 and 0.29 x 10 22 cm 2 . In the first case, the first temperature is fixed to 0.617 keV, whereas in the 
second case, both temperatures are fixed ( cooler component at 0.617 and hotter component at 1.89 keV). Fluxes refer to the 0.5-10.0 keV band, except for ROSAT data, where they are given in the 
0.5-2.0 keV band. 
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Both PSPC datasets were taken far from periastron, at phases 
0.41 (in 1991) and 0.28 (in 1993), so that we expected few 
changes in recorded intensity. Indeed, the count rates, fluxes, 
and spectral parameters of the two observations agree with each 
other within the errors. The 1993 data was also taken at a similar 
phase as the sixth XMM-Newton observation, and there is also 
fair agreement between them. HRI observations occurred just 
before and just after a periastron passage at phases 0.95 for the 

1994 observation (i.e., similar in phase to the third Swift obser- 
vation) and 0.11 for the 1995 observation (i.e., similar in phase 
to the fifth Swift observation). In fact, the recorded count rates 
are higher in the 1994 data, which is closest to the periastron 
time. 

Using the HRI response matrix and the same model as in the 
last XMM-Newton observation (i.e., same temperatures, absorp- 
tions, and normalization ratio), we found an observed flux in the 
0.5-2.keV band of 8 and 5 x 10" 13 erg cirT 2 s _1 for the 1994 and 

1995 observations, respectively. These fluxes agree well with the 
ones observed at similar phases with Swift. Overall, the ROSAT 
data thus appear compatible with the variations seen in XMM- 
Newton and Swift observations. 

ASCA observed Cyg OB2 in April 1993, at phase 0.27, nearly 
at the same time as the second ROSAT-PSPC observation. We 
extracted SIS spectra for both SIS-0 and SIS-1, and for both 
bright and bright2 modes (i.e., 4 spectra in total, implying that 
we merged data taken with different bit rates). For the source, 
we used the same region as in Swift and ROSAT, whereas the 
background region was chosen as a nearby circle of the same ra- 
dius and on the same CCD. The source region is quite small for 
the ASCA resolution, but Cyg OB2 #9 appears close to a CCD 
edge. Another problem in this observation is the high back- 
ground: first, straylight from the nearby extremely bright source 
Cyg X-3 is contaminating the field-of-view; second, Cyg OB2 #9 
appears on the PSF wings of the bright Cyg OB2 #8 A. Indeed, 
Cyg OB2 #9 is not detectable by eye in single datasets, and it ap- 
pears only as a faint object in the SIS global image. The result- 
ing spectra are thus noisy, but we nevertheless tried to fit them 
with the same model as above (with temperatures and absorption 
fixed, see Tabled. The observed flux, in the 0.8-4. keV range 
amounts to (8.4 + 0.5) x 10~ 13 ergcrcT 2 s _1 , slightly smaller but 
within the e rrors of the value of ' ( 1.1 + 0.4) x 10 _12 ergcm _2 s _1 
reported by Kitam oto & Mukail (Q996). Compared to the sixth 
XMM-Newton dataset (obtained at a similar phase) or the con- 
temporary ROSAT observation, we again find a smaller flux in 
our ASCA fitting. However, considering the background and res- 
olution problems, the agreement is not that bad. 

In summary, the archival data agree well with the XMM- 
Newton and Swift results, providing further support to the pres- 
ence of phase-locked variations in Cyg OB2 #9 since these dif- 
ferent datasets were taken several periods apart. 

4.3. A wind-wind collision in CygOB2#9 

The presence of bright and hard X-ray emission (with kT > 
1.5 keV) is quite unusual in O-type stars. It may be linked to 
a magnetically-confined wind in a single star or to a wind- 
wind collision in a binary system. In our case, even though the 
nonthermal radio emission implies the presence of a magnetic 
field, the latter explanation seems to be the correct one. Indeed, 
Cyg OB2 #9 displays long-term variations, apparently compati- 
ble with the timescale of the orbital period. 

Theoretical considerations show that there are two main 
sources of phase-locked va riability am ongst colliding win d bi- 
naries dStevens etalJll992l see also Giidel & Nazl l2009l for 



Table 5. Distance and flux ratios for XMM-Newton and Swift 
data, for the fits with the second temperature free to vary and 
with the observations taken in Oct04 and in Janl 1 as references, 
respectively. 



Date 


D ref /D 




F X /F; 


W 










0.5-10. keV 


0.5-2. keV 


2.-10. keV 


XMM-Newton 












Nov 04a 


0.96 


0.93 ± 0.03 


0.93 ± 


0.03 


0.95 


±0.04 


Nov 04b 


0.93 


1.04 ±0.03 


1.06 ± 


0.03 


0.99 


±0.04 


Nov 04c 


0.90 


0.95 ± 0.04 


0.92 ± 


0.04 


1.01 


±0.06 


Apr 07 


0.84 


0.83 ± 0.03 


0.80 ± 


0.03 


0.91 


±0.04 


May 07 


0.83 


0.84 ± 0.03 


0.82 ± 


0.03 


0.91 


±0.05 


Jun 11 


3.90 


3.97 ± 0.09 


4.17 ± 


0.09 


3.43 


±0.11 


Swift 














Apr 11 


1.51 


1.32 ±0.19 


1.21 ± 


0.22 


1.56 


±0.32 


May 11 


2.43 


3.14 ±0.37 


3.39 ± 


0.49 


2.62 


±0.42 


Jul 11 


4.27 


4.01 ±0.51 


4.46 ± 


0.68 


3.02 


±0.55 


Oct 11 


1.35 


1.90 ±0.30 


2.11 ± 


0.36 


1.48 


±0.29 



a review of recent observations). The first one is the absorp- 
tion. Since the collision is seen alternatively through the wind 
of each component, changes in absorption are expected. This is 
especially relevant for systems with two very different winds, 
e.g. WR+O: two spectacular cases in this respe ct are y 1 Vel 
(Will is et alll995l) and WR22 dGosset et al.ll2009l) . The similar- 
ity of the two components in Cyg OB2 #9 makes such variations 
less prominent; indeed, no significant change in absorption is 
detected in the data. 

The second source of phase-lo cked variability is th e chang- 
ing separation in eccentric binaries dStevens et al.ll992b . For adi- 
abatic cases, the separation between the stars is high, and the 
X-ray luminosity is expected to follow a 1/D relation (i.e., it 
should be maximum at periastron) with a minimal impact from 
the velocity changes because the winds are close to their maxi- 
mal speeds at all phases. For radiative cases, the X-ray luminos- 
ity and plasma temperatures should be minimum at periastron, as 
the winds then collide at lower speeds. Adiabatic systems should 
be those with long periods, like Cyg OB2 #9, but not many long- 
period binaries have been monitored in X-rays, and results were 
mixed, y 1 Vel and WR22 did not show any 1 j D variation, despite 
the changing separ ation of their components (Rauw et"aTI 2000b; 
iGosset et alj|2009l) . CygOB2#8A and WR140 deviate from ex- 
pecta tions at periastron, since the collisions then become radia- 
tive dDe Becke r et al . 20061 iPollock et al. 20021 iDe Becke r et alJ 
|2012|) . Hints of the influence of a 1/D relation may h ave been 
found in HD 93205 (see Fig. 2 in lAntokhin et alj|2003l) . though 
the results need confirmation, in view of the huge noise, and 
in HD 93403, though the intrinsic emission of the stars blurred 
the results, yielding a luminosity change lower than that ex- 
pected from a simple 1/D relation dRauw et al.ll2002l) . Among 
evolved systems, WR25 se ems to display a rather "clean" 1 /D 
variation in its luminosity dGosset 11200 7; Giidel & Naze 2009). 
For m ost of its orbit, it is also the case of rj Car (Okazaki et al. 
2008), though the system shows large cycle-to-cycle variations 
in flux, and there are strong indications that the wind collision re- 
gion co llapse onto the companion star at periastron ( Parki n et all 
1201 Id and references therein). Up to now, no comparable exam- 
ple among O+O binaries has been found, so Cyg OB2 #9 may 
well become the first one. 

Indeed, while there seems to be no change in absorption, 
at least for our sampling, the normalization factors and fluxes 
clearly change with time. They show a sharp increase near peri- 
astron (see middle panel of Fig. [8). In fact, a strong correlation 
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Fig. 8. Variations in the spectral parameters (left) and observed flux (middle) with phase, for the best-fit model with the first tem- 
perature and the absorptions fixed. The solid and open symbols correspond to XMM-Newton and Swift data, respectively. The right 
panel shows, for the same model, the flux ratio F^ abs / Ff^j! as a function of the inverse stellar separation ratio D re f/D for the 
XMM-Newton and Swift data. The reference data corresponds to the first observation of each dataset, i.e., data from Oct 04 for 
XMM-Newton and data from Jan 1 1 for Swift. The black circles corresponds to the total, 0.5-10.0keV flux, the green squares to the 
soft, 0.5-2.0keV flux, and the blue triangles to the hard, 2.0-10.0keV flux. 



appears when the fluxes are compared to the inverse of the stel- 
lar separation. The correlation is not perfect, but we know that 
there are instrumental differences (up to 20% in flux, see above) 
and larger errors for the less sensitive Swift data. To reduce the 
impact of these instrumental effects, we performed a comparison 
within each dataset (see last panel of Fig.[8]and Table|5}: within 
the errors, the fluxes vary as 1 I D, as expected from an adiabatic 
collision in an eccentric system dStevens et al.ll 1 992b . In fact, the 
linear correlation coefficient between the distance and flux ratios 
listed in Tablets 0.97 for Swift +XMM-Newton data and 0.999 
for the more precise XMM-Newton data alone. In this context, 
it may also be noted that the increase by a factor of 1 .6 in the 
ROSAT-HR1 count rates is well matched by the decrease in sep- 
aration between orbital phase = 0.11 and <p — 0.95 by the same 
factor. The only slight departure from the 1 /D relation (Fig. [8) is 
the hard X-ray flux appearing slightly below the expected value: 
this may be due to the collision becoming slightly radiative be- 
cause of the increased density at periastron. This would also ex- 
plain the presence of cool plasma responsible of Ha absorption 
at these phases in the optical spectra. However, this does not af- 
fect the bulk of the X-ray emission, since the l/D correlation is 
closely followed overall, with no evidence of blurring by an ad- 
ditional, constant intrinsic emission from the stars. The absence 
of dramatic changes in Ha also correlates well with an adiabatic 
nature for the wind-wind collision. 

On the other hand, the second temperature of the spectral 
fitting is clearly lower near periastron in the XMM-Newton 
data (see left panel of Fig. [8}. This may be explained by a 
shock occurring at a smaller separation. Indeed, a shorter dis- 
tance between the two stars may imply that winds collide be- 
fore they could accelerate to full speed, resulting in a cooler 
shocked plasm a. Moreover, at small separations, radiative brak- 
ing/inhibition dStevens & Pollockll994tlGavlev et all 19971) may 
take place, slowing down the winds further before they collide. 
Such a cooling of the hard component is predi cted by models 
of win d-wind collisions in eccentric binaries: iPittard & Parkin! 
d2010l see their Fig. 22) find a change by a factor of 3 in the 
highest temperature, whereas Parkin & Gossetl (1201 lal) find an 



order of magnitude variation when the shock plunges deep into 
the wind acceleration region. In CygOB2#9, the temperature 
only changes by 25%, but the separatio n is m uch larger than 
in the case modeled by Pittard & Parkin (2010) and the wind- 
wi nd collision is closer to balance than in the case modeled 
bv lParkin & Gossetl d201 lal) . Therefore, we expect the influence 
of wind acceleration in Cyg OB2 #9 to be less pronounced than 
found by these authors, weakening the cooling effects. 

To put these observing results in perspective, we did a first 
modeling of the wind- wind collision in CygOB 2#9. We here 
use the analytical models of lStevens et al.l (1 19921 see als o a de- 
tailed example of the calculation in iParkin et al.ll201 lbl) . From 
the ram pressure balance, they allow us to find the position of the 
collision, the cone's opening angle, plasma temperatures (from 
the pre-shock wind velocitiefl), and the cooling parameter (to 
determine the nature of the collision, adiabatic if > 1, radia- 
tive otherwise). As input, such models require wind and stel- 
lar parameters: stellar radii and separation, wind mass-loss rates, 
and velocities. The orbital solution has just been found, yielding 
the separation at each phase. However, not much information is 
available on the individual wind parameters, only approximate 
spectral types have been derived. We therefore made the follow- 
ing assumptions: inclination of the system of about 62 degrees 
(so that the true ma sses are compatible with the predictions of 
Marti ns et al.ll2005l) . typi cal stellar parameters of 05-5.51 and 
03-4III stars taken from Martin s et al.l d2005l). a nd mass-loss 
rates derived from the recipe of IVink et al.l d2000t) . This yields 
mass-loss rates of 5 x 10~ 6 M yr" 1 and 1.2 x 10~ 5 M yr _1 and 
terminal velocities (= 2.6xescape velocities) of 2540 km s~' and 
2440km s _1 for the primary and secondary, respectively. We also 
assume the usual /3-law with fj = 1 for the wind velocity evolu- 
tion. At maximum separation, this results in a wind momentum 
ratio of 0.43 (or 2.3 in favor of the secondary), which gives an 



9 We assume that the emission-weighted post-shock temperature is 
kT(in keV) = 0.6 x v(in 1000 kms" 1 ) 2 . This is roughly half the tem- 
perature at the stagnation point and is meant to account for the lower 
temperatures obtained in the downstream flow as the off-axis shocks 
become oblique. 
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opening angle of the collision zone of abo ut 150 degrees, com- 
patibl e with the radio map (see Fig. 6 of iDoughertv & Pittard! 
[20Q6L note that the angle estimated from that figure is approxi- 
mate, because sensivity or wind opacity may affect its determi- 
nation). Figure [9] shows the evolution with phase of several pa- 
rameters of the collision zone. The position of the balance point 
(at 0.4D from the pri mary) does not chang e with phase, and the 
cooling parameter ;f dStevens et al.lll992h is always >> 1, im- 
plying that the shocks should still be relatively adiabatic even 
at periastron (which nicely conforms with the X-ray observa- 
tions). The agreement is not perfect, however, since the shocked 
plasma temperatures appear too high (about 3.6 keV at apas- 
tron and 3.1keV at periastron vs observed values of 2.5 keV 
and 1 .9 keV), and the relative change in temperatures is slightly 
too low (13% vs 25% observed). We therefore include the ef- 
fect of radiative inhibition in the model: rather than using a sim- 
ple /3 velocity law, we here assume the winds to be radiatively 
driven following the ICastor et alj d 1975b approximation and in- 
cluding the influence of the opposing star radiation field (e.g. 
IStevens & Poilocl[T99l: iGavlev et alJll997h . Such a model im- 
proves the agreement with observations: the winds collide at 
lower speeds, even at apastron, resulting in lower plasma tem- 
peratures (about 2.7 keV at apastron and 2.0 keV at periastron), 
and the decrease in plasma temperature is now 25% between 
apastron and periastron (Fig. [§)■ In this case, the shock still re- 
mains mostly adiabatic throughout the orbit, as observed. Wind 
clum ping may play a further role in decreasing plasma tempera- 
tures (Pittard 2007, and references therein). The next step will be 
to make full 3D hydrodynamical simulations of th e wind-wind 
collision (c.f.Hlttard 2009: lParkin& Gosset 201 lab , to calculate 
synthetic spectra and lightcurves. Their comparison with obser- 
vations will allow determination of the density, temperature, and 
velocity of the hot plasma within the wind-wind collision region, 
helping to constrain the stellar mass-loss rates and wind speeds 
in this system. Constraints on the magnetic field within the colli- 
sion region will also be derived. However, this detailed modeling 
requires all data to be at hand (not only optical and X-ray obser- 
vations - see future papers of the series); it is also clearly beyond 
the scope of this paper. 

5. Summary and conclusions 

In mid-2011 the first observable periastron passage of 
CygOB2#9 occurred since the binarity discovery. To observe 
this event under all aspects, we organized a multiwavelength 
campaign with a simultaneous coverage at radio, infrared, op- 
tical, and X-ray energies. Here we report on the results of the 
optical and X-ray monitorings. 

The optical spectroscopic monitoring obtained with the 
Sophie spectrograph enabled us to refine the orbital solution, 
confirming preliminary results. A spectral disentangling was 
also performed, revealing the two stars to be of 03-4III and 05- 
5.51 types for the secondary and primary, respectively. While the 
He ii 14686 A line follows the motion of the hottest component 
(though with a shift), the Ha line appears extremely broad but 
unshifted. The line profile changes with phase, with both en- 
hanced absorption and stronger emission at periastron. 

X-ray observations by XMM-Newton and Swift have been 
obtained as Cyg OB2 #9 reached the 201 1 periastron. The high- 
energy monitoring revealed a rather constant X-ray spectrum, 
except at periastron where the emission becomes slightly softer. 
The second thermal component then has a temperature reduced 
by about 25% at periastron, which is reproduced well by colli- 
sion models including radiative inhibition. The monitoring also 
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Fig. 9. Evolution with phase of the collision parameters (position 
of balance point, cooling parameter^, pre-shock velocities, and 
average plasma temperature, kT - 0.6 x v(in 1000 km s -1 ) 2 ). 
Thick red lines and thin black lines correspond to the cases with 
and without radiative inhibition, respectively. 



unveiled the luminosity variations of the system, which follow 
the l/D changes predicted for adiabatic collisions remarkably 
closely. The only (slight) departure from that relation is a very 
small decrease in the hard-band flux, which could indicate that 
the collision has become very slightly radiative at periastron. 
Overall, CygOB2#9 seems the best (and the first!) example of 
an adiabatic O+O colliding-wind system. 
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Table 1. Journal of observations. 
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Fig. 1. The spectrum from June 2010 (black thin line) and the 
DIB correction (red thick line). 



Date 


HJD-2 450 000 





Sophie 






Sep 07 


4348.765 


0.379 


Oct 07 


4379.063 


0.414 


Dec 07 


4463.748 


0.512 


Mar 09 


4906.675 


0.028 


Apr 09 


4936.968 


0.064 


May 09 


4953.595 


0.083 


Jun 09 


5006.497 


0.145 


Jul 09 


5022.578 


0.164 


Aug 09 


5058.433 


0.205 


Sep 09 


5079.965 


0.230 


Oct 09 


5123.370 


0.281 


Nov 09 


5142.240 


0.303 


Apr 10 


5316.621 


0.506 


May 10 


5346.054 


0.540 


Jun 10 


5372.560 


0.571 


Jul 10 


5400.557 


0.604 


Aug 10 


5433.963 


0.643 


Sep 10 


5444.420 


0.655 


Oct 10 


5482.870 


0.700 


Nov 10 


5523.761 


0.747 


Dec 10 


5557.230 


0.786 


Jan 11 


5584.273 


0.818 


Mar 11 


5626.701 


0.867 


Apr 1 1 


5675.632 


0.924 


May 11 


5698.603 


0.951 


Jun 11a 


5731.576 


0.989 


Jun lib 


5737.249 


0.996 


Jul 11a 


5745.042 


0.005 


Jul 1 lb 


5755.550 


0.017 


Jul 1 lc 


5769.524 


0.034 


Aug 11 


5790.465 


0.058 


WIRO (2011) 






May 31 


5713.847 


0.969 


Jun 01 


5714.763 


0.970 


Jun 03 


5716.737 


0.972 


Jun 04 


5717.742 


0.973 


Jun 05 


5718.738 


0.975 


Jun 13 


5726.765 


0.984 


Jun 14 


5727.794 


0.985 


Jun 24 


5737.758 


0.997 


Jun 25 


5738.778 


0.998 


Jun 26 


5739.444 


0.999 


Jun 27 


5740.441 


0.000 


Jun 28 


5741.443 


0.001 


Jul 15 


5758.752 


0.021 


Jul 23 


5766.924 


0.031 


Jul 24 


5767.438 


0.031 


Jul 25 


5768.452 


0.032 


XMM-Newton 






Oct 04 


T)AO con 

3308.580 


0.167 


Nov 04a 


3318.558 


0.178 


Nov 04b 


3328.544 


0.190 


Nov 04c 


3338.506 


0.202 


Apr 07 


4220.355 


0.229 


A/fot; (Y7 

iviay u/ 


4ZZ4. 1 / U 


n oil 


Jun 11 


5738.255 


0.997 


Swift 






Jan 11 


5571.618 


0.803 


Apr 11 


5655.837 


0.901 


May 11 


5700.082 


0.953 


Jul 11 


5743.839 


0.004 


Oct 11 


5842.169 


0.118 



Notes. Datasets are identified by the month of observation followed by 
the last two digits of the year, except for the WIRO data where the 
numbers correspond to the dates. 
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Fig. 4. RVs determined by simple fits (top) and by finding the best shift of fixed Gaussians (bottom left). Small cyan symbols 
correspond to WIRO data, other symbols are the same in all panels, with one colored symbol for each line. The RV measured for 
a narrow interstellar feature is superimposed on the Hen 14686 A subpanel using asterisks, to show the quality of the wavelength 
calibration. The bottom right panel shows the RVs determined by disentangling. The best-fit orbital solution (Table [3} found from 
disentangling is superimposed on each panel. 
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Fig. 5. Disentangling results. Each panel shows, from top to bottom, the spectra of the primary star, the secondary star, and the 
interstellar features (with continuum shifted to +1) in four different spectral domains. The normalization is imperfect, as is often the 
case for disentangling results when the orbital RV shifts are modest. 
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Table 2. Radial velocities found from disentangling for the 
Sophie data. 



Date 


HJD-2450 000 


RVi 


RV 2 






km s -1 


km s -1 


Sep 07 


4348.765 


-14.2 


-11.1 


Oct 07 


4379.063 


-18.3 


-17.9 


Dec 07 


4463.748 


-14.5 


-9.3 


Mar 09 


4906.675 


-83.9 


60.4 


Apr 09 


4936.968 


-32.1 


24.3 


May 09 


4953.595 


-38.1 


-7.9 


Jun 09 


5006.497 


-47.9 


-11.6 


Jul 09 


5022.578 


-28.9 


-13.5 


Aug 09 


5058.433 


-15.7 


-16.0 


Sep 09 


5079.965 


-24.6 


-21.4 


Oct 09 


5123.370 


-12.6 


-26.0 


Nov 09 


5142.240 


-17.8 


-16.9 


Apr 10 


5316.621 


-4.0 


-33.0 


May 10 


5346.054 


-9.5 


-27.4 


Jun 10 


5372.560 


-17.4 


-15.8 


Jul 10 


5400.557 


-24.2 


-20.5 


Aug 10 


5433.963 


-23.9 


-19.9 


Sep 10 


5444.420 


-16.2 


-12.5 


Oct 10 


5482.870 


-15.1 


-9.6 


Nov 10 


5523.761 


-38.9 


-6.8 


Dec 10 


5557.230 


-37.8 


2.0 


Jan 11 


5584.273 


-48.2 


-5.9 


Mar 11 


5626.701 


-51.4 


28.1 


Apr 11 


5675.632 


-53.9 


41.2 


May 11 


5698.603 


-68.7 


55.8 


Jun 11a 


5731.576 


-140.9 


117.4 


Jun lib 


5737.249 


-143.8 


112.5 


Jul 11a 


5745.042 


-140.9 


96.7 


Jul lib 


5755.550 


-90.3 


95.6 


Jul 11c 


5769.524 


-88.0 


48.2 


Aug 11 


5790.465 


-32.2 


39.1 



